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We study by Monte Carlo simulation a binary mixture of neutral and dipolar hard- 
spheres with non-additive diameters. With a view to understanding the interplay 
between population inversion for an open pore and the demixing phase transitions, 
the mixture is considered in the bulk and confined between two parallel hard-walls 
modeling a slit pore. A uniform field is applied in the pore in order to control its 
composition as shown previously. The demixing lines in the bulk and in the pore are 
studied by the Gibbs Ensemble Monte Carlo method. The open pore-bulk mixture 
equilibrium is studied by a combination of canonical/grand canonical simulations. 
A moderate electrostatic coupling is considered for remaining close to the conditions 
in which a jump in the adsorption of the minority species has been observed at 
zero field. Demixing lines are given in the bulk and for two different pore widths in 
parallel and normal fields, together with population inversion paths. Similarly to the 
effect of geometrical confinement, a normal field is found to favor the mixed state so 
that the population inversion does not interfere with demixing. A parallel field leads 
to more complex scenarios. Some indications for future work are briefly discussed. 

I. INTRODUCTION 

Fluid mixtures adsorbed in porous materials have been studied in the past decades with 
motivations ranging from technological problems to answering basic theoretical questions 
(see for e.g. [ll, Q])- Understanding the factors that regulate the behavior of the adsorbed 
fluid such as its composition, the adsorption geometry, the different interactions, etc. is 
thus of central importance. Among the various mixtures that are studied, colloidal ones 
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play a special role since it is possible to tune to some extent (eg. by adding polymeric 
depletants) the effective interaction between the particles and their interaction with the 
confining medium js, ^. Since colloidal particles are also often polar or polarizable, the 
combination of geometrical confinement and the action of an external field offers then a 
unique way to control their behavior in confinement. In recent work, we have shown [sl, 
this by Monte Carlo simulation for the simple model of a binary mixture of neutral and 
dipolar hard-spheres confined between two parallel hard-walls modeling a slit pore. We 
have shown that the composition of the adsorbed fluid can be controlled in a flexible way 
by applying an external field in the pore, at fixed thermodynamic state of the bulk fluid. 
While the effect of an external field on confined polar fluids is well known (see for example 

n fl 

references in [7|] for theoretical work and [8| for experiments) the originality of the method 
we have proposed to control the composition is the combination of the demixing instability 
in a mixture with unfavorable attraction (or greater repulsion) between unlike particles and 
the lowering of the potential energy of the dipolar fluid in the region where the field is 
applied. This also differs from the numerous studies in which the adsorption is controlled 
through some action on the bulk fluid such a change in pressure. The possibility to control 
the behavior of the confined fluid at fixed bulk state would indeed be a great advantage 
for some applications. Not having to rely on the subtle combination of specific interactions 
should also result in a robust method, feasible with simple components. This is why we 
choose to illustrate the basic mechanisms js, 6| on the simple mixture of neutral and dipolar 
hard-spheres and perfectly smooth planar walls. The physical situation closest to this model 
is then the adsorption of colloids, in contrast with molecular adsorption in which specific 
interactions should be considered. The dipolar hard-spheres and the neutral ones interact 
with a hard-sphere potential with non-additive diameter. The excess repulsion associated 
with a positive non-additivity is known to favor demixing {9]. This is one of the basic 
mechanisms for observing the population inversion in the pore-bulk equilibrium. The main 
result we have obtained in refs. is the possibility to produce, in the pore, a field 

induced population inversion when the bulk is near the demixing instability (referred to as 
the f-PINBI effect): we start from the situation of a pore in equilibrium with a bulk mixture 
having a tendency to demix and widely different concentrations of the two components; 
when the strength of the field applied in the pore region is varied, a jump in the adsorption 
of the minority component occurs. After a certain threshold value, the pore becomes filled 
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with this species, while the former majority component simultaneously desorbs. One may 
accordingly produce a jump in the physical properties related to the fluid cornposition by 



the sole action of the external fleld. A similar effect has been described in refs lHj , but 
it was induced by a variation of the bulk density. Here the thermodynamic state of the bulk 
fluid is unaffected. 

Both for practical and theoretical reasons, it is important to understand the connection 
between this phenomenon of population inversion and the phase transitions in the bulk and 
in the conflned fluid. For speciflc applications, one may indeed wish to know in advance 
how to tune the controllable parameters in order to favor the population inversion with 
respect to the demixing in the pore and vice versa. More fundamentally, the phase diagram 
of binary mixtures is known to be quite complex, even in the simplest case with symmetric 



interactions - see for instance 



12| for a general discussion and 
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4l5l | for explicit calculations 



on model systems showing for example the progressive substitution of the evaporation- 
condensation by demixing as temperature increases. With dipolar species, the presence 
of additional parameters (dipole moment, fleld strength) complicates further the general 
phase diagram of the mixture, already in the bulk - see for instance refs 16l-l21|. The 
conflnement increases the diversity of possible scenarios, including speciflc effects such as 
layering, capillary condensation, etc - see for e.g. 22N24| . This domain is accordingly even 



less known than the phase transitions of bulk mixtures. The theoretical study by semi- 
analytic methods of such inhomogeneous fluids with anisotropic interactions is accordingly 
also problematic (see 25| for a recent study of conflned dipoles). 

The objective of this work is thus to analyze further this question on the same model as in 
ref. {5,6]. We will supplement our previous results relative to the bulk phase diagram of the 
non-additive hard-sphere mixture and the population inversion in the presence of dipoles by 
the phase separation in the bulk and in the pore in the presence of dipoles, flrst at zero fleld. 
Field strength-composition coexistence curves in the pore will then be given for selected 
values of the total density, with no attempt to do a full mapping of the parameters space. 
Indeed, even in this simple model, we have at least four physical parameters: non-additivity 
parameter 6, pore width, fleld strength (parallel or normal), dipole moment or temperature, 
besides the usual composition variables. It is thus clear than only a limited domain of the 
parameters space can be explored by simulation in a reasonable time. For the state points 
considered, we will try to clarify the status of the population inversion with respect to the 
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demixing transitions in the bulk and in the pore. 

Another question is the choice of the control variables since this determines the simulation 
strategy for the bulk-bulk, pore-pore and pore-bulk equilibria. The first one was determined 
by the standard Gibbs Ensemble Monte Carlo method of Panagiotopoulos 26[ (GEMC). 
The pore-pore equilibrium was studied by the appropriate version of this method. Even 
if the open pore for which the PINBI effect occurs can be studied similarly, we wish to use as 
the control variables the total density and mole fractions in the bulk. These control variables 
are indeed those that are easily accessible experimentally. We thus adopted a strategy which 
combines a canonical/grand canonical simulation: the chemical potentials are measured in 
the bulk considered in the canonical ensemble. They are next used to study the fluid in the 
open pore in the grand canonical (GC) ensemble. 

This paper is organized as follows: In section II, we present the model and detail the 
simulation methods we used to study it in the different situation (bulk, pore and bulk- 
pore equilibria). In section III, we present some corresponding coexistence curves, without 
field. We then discuss the relation between the field-induced population inversion and the 
coexistence for selected values of the total density. In the conclusion, we summarize the 
main results and indicate some directions for future work. 



II. MODELS AND SIMULATION TECHNIQUES 



A. Model 



The model considered is a binary mixture of neutral hard-spheres (label 1) and dipolar 
hard-spheres (label 2) of equal diameter a, with a non-additive diameter ai2 = cr(l + 6). 
This parameter is kept fixed to 6 = 0.2 in this study. The dipolar hard-spheres bare a point 
dipole fi at their center. In addition to the hard-sphere interaction, the dipoles interact then 
through the dipolar potential: 

^ fMfJ^ _ g (^i.rij)(^j.rij) 

Dip »,3 j,5 ■ V / 

ij ij 

The reduced dipole moment is given by /i* = n/ ^/ksTa^ . For this study, we took /i* = 1. 
With this magnitude, an energy ^*E* = 1 would correspond for colloids with diameter 
1/im, for example, to /i = 2 lO^D and E* = 1 to E ^ 6.1 10~^V/fim at T = 300K. 
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This relatively weak value has been chosen so as to remain close to the situation considered 
in our previous work s], O]: when a non-additive mixture of neutral hard-spheres with 
6 = 0.2 is confined in a pore of width H = 3ct a density induced PINBI occurs for a 
bulk state {p^ ^ 0.555, X2 = 0.02). This non-additivity of the hard-sphere interaction is a 
distinguishing feature of the model considered here with respect to the mixture of neutral 
and dipolar hard-spheres considered in previous studies. 

For the purpose of understanding the interplay between the properties of the fiuid in the 
bulk and in confinement - here the connection between population inversion for an open 
pore and the demixing phase transitions in the bulk and in the confined fiuid - this model 
will be studied in three steps for which we describe below the simulation method used: the 
bulk-bulk, the pore-pore and bulk-pore equilibria. 

B. Simulation of the bulk fluid 

To study the mixing/demixing transition for the bulk fiuid, we performed standard Gibbs 



MC simulation in the (A'^i, ^t? T) ensemble |26|. To achieve this, one has to simulate two 
cubic boxes {A and B) of fixed total volume Vr = Va + Vb with periodic boundary conditions 
in each direction. The boxes exchange particles and volume (the box sizes La and Lb change 



during the simulation 



sums with conducting boundary conditions 
took the same value as in our previous work 



26|). The dipolar energy of the periodic system is evaluated by Ewald 

. For the inverse screening parameter a we 
6] a = 7 / La,b The sum in reciprocal space 
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included all lattice vectors k = n = {nx^ny^riz), with < 80. 

The control variables are the reduced total density = {Ni + N2)<j'^ /Vt and temperature 
defined as T* = kept fixed here to T* = 1. To study the coexistence at a given p\, 

we started the simulations with box lengths ranging from La = Lb = 12cr for nearly mixed 
fiuids to La = Lb = 15.5cr. This keeps a representative number of the minority species 
when the system is nearly completely demixed. In all cases, we took Ni = N2 particles, a 
choice appropriate to the moderately asymmetric coexistence curves considered here. Each 
MC cycle consisted in 5A^t' particle exchanges, Nj- particle translation moves, N2 rotation 
moves and 1 volume exchange. The type of movement was chosen at random. We performed 
typically 10000 equilibration cycles and from 2^^ to 2^^ for critical conditions of accumulation 
cycles. Note that because of the weak dipolar interaction, dipole exchanges remain efficient. 
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For lower temperature or higher dipole moment, dipole exchanges become difficult (in this 



17( 1 but this was not 



case one might use the identity change proposed by Blair and Patey 

necessary in the conditions we have investigated). 

When the mixture separates, one measures the mean total density in each box: p\ =< 
^ — > and its mean composition in species i : =< j^a^^a > (and same for box 

B). Note that near the critical demixing point the phases exchange between the two boxes. 

The equilibrium is then determined from the compositions histograms. 



C. Simulation of the confined fiuid 

The fluid is now confined between two parallel hard walls modeling a slit pore of width H. 
Due to the hard-core interaction with the walls, the centers of the particles of both species 
can move within the interval [a/2, H — a/2]. An electric field either normal E = E or 
parallel E = E Uj. to the walls is applied in the pore. To avoid the complication due to 
image forces, we assume continuity of the dielectric constant at the wall boundary. Since this 
excludes the case of conducting walls, the field is assumed to be generated by some external 
device, say parallel plates well away from the walls. To take into account the confinement 
one has to modify the Ewald summation. Periodic boundary conditions are now applied 
only in the x and y directions. We used the simple modification proposed by Yeh and 
Berkovitz 



29|: a simulation box of dimensions L^, Ly, is used in which a vacuum region 
is incorporated in the z-direction such a.s Lz = H + /^vacuum = iL where L is the lateral 
dimension: = Ly = L. Our implementation of this modified 3D-Ewald summation is 



similar to that detailed in the paper by Klapp and Schoen [30|. We took for the Ewald 
summation the set of parameters: a = 7/L where L is the lateral box size, 7 = 10 that gives 
a gap width 7L — H and = 80, with n = {ux, Uy, riz/'y), for the terms in reciprocal space. 



A large gap might be avoided by using the electrostatic layer correction (ELC) 3l|, |32|. We 
did not try this here since our objective was to test the physical ideas rather than to propose 
a fully optimized algorithm. 

To ensure the pore-pore equilibrium, the GEMC simulation must be conducted as shown 
in the volume exchange being replaced by surface exchange. Simulations are thus 
performed with two slits of fixed width H, fixed total surface Sa + Sb , with Sa,b = L\ ^ 
and fixed total number of particles. In the corresponding {Ni,N2,Vt,T) Gibbs ensemble 
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Vt = {Sa + Sb){H — a) is the volume accessible to the particles centers, that is a/2 away 
from the walls. We thus kept the same control parameters as for the coexistence in the bulk 
: = (A''i + N2)(t'^/Vt (the superscript p being used only to distinguish the pore from the 
bulk). However, we do not impose the equilibrium with the bulk fluid. The first reason is 
a practical one since the GEMC simulation does not require the knowledge of the chemical 
potentials. The second one is that a closed and an open pore are two physically different 
situations. In particular, equal densities in the pore and in the bulk does not imply the 
equality of the chemical potentials. For an open pore, the average density is imposed by the 
equality of the chemical potentials through the exchange of particles with the bulk. If that 
case, the correspondence between in the bulk and the average density p in the pore can 
be done by determining the chemical potentials from a (A^, V, T) simulation in the bulk and 
a GC one in the pore. Finally, the reduced temperature is kept fixed to T* = 1. Isothermal 
paths at varying field strength correspond then to a fixed /i*. 

As in the bulk, the mean total density in box A for the phase separated confined fluid is 
p\ =< ''^^ ~^v^ ^ > s-iid its mean composition in species i is =< j^a^j^a > (and same 
for box B). 

For H = 3(7, we start the simulations with boxes of equal sizes ranging from = Lb = 
20a to La = Lb = 21.5a for nearly completely demixed fluids. For H = 9a, we start 
with La = Lb = I2a and La = Lb = 15cr for highly asymmetric composition. Besides 
these points that are specific to the confinement geometry, the dynamics of the simulation 
(number of cycles, etc.) was the same as in the bulk studies. It should be noticed that 
fluctuations of the volume and of the particle numbers can be significant during the Gibbs 
ensemble simulation, due to the limitation of the number of dipoles that can be taken for 
a reasonable restitution time. It is thus difficult to accurately determine the coexistence 
points (see the irregularities in the slopes of the tie lines ) but we have checked that longer 
runs did not significantly change the coexistence lines. 

Finally it should be noticed that such finite size effects become significant when ap- 
proaching the critical point of the mixing/ demixing transition. As an example, a symmetric 
mixture of non-additive hard-spheres with p^ = 0.46 confined in a pore of width H = 9a is 
found to be slightly demixed (with = 0.46) for a box size L = 12a, whereas it remains 
mixed (with = 0.5) for L = 18a. It is difficult to perform such studies with dipolar hard- 
spheres due to the number of dipoles to be considered when the boxes are large, especially 
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in confinement. Since locating precisely the critical point is not our goal here, we did not 
investigate this point further. 



D. Fluid confined in an open pore 

To establish the connection between the population inversion in a mixture confined in 
an open pore and the demixing lines in the bulk and in the pore, we recall first the method 
used in our previous work to obtain the population inversion points. The pore-bulk mixture 
equihbrium is determined by the equality of the chemical potentials fii and /i2, in the 



bulk and in the pore. As in ref. jsl we used for this a combined canonical/grand canonical 
simulation. Using a MC simulation in the (A'^i, A^2, V, T) canonical ensemble, we first measure 



by Widom's insertion method 
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the chemical potential of both species. The bulk fluid 
state is then characterized by a total density and composition Xg. Note that these 
overall density and composition need not necessarily correspond to a homogeneous mixture. 
The chemical potentials /ii,2(PT' ^2) obtained are next imposed to the slab through a GC 
simulation from which one measures the mean density p{p^, Xg) and the mean composition in 
the pore X2{p^, Xg). In the density-induced PINBI, the population inversion is characterized 
by a set of X2(px) curves, one for each value of X2 considered. They describe the change of 
the mean pore composition from the bulk one, For the f-PINBI instead, the inversion 
with respect to Xg is characterized by a set of X2{E*) curves, all for the same p^. 



III. RESULTS, DISCUSSION 
A. Bulk phase diagram 

Fig. [1] compares the bulk phase diagram of a mixture of neutral and dipolar hard- 
spheres with that of the hard-spheres mixture without dipoles. In both cases, the non- 
additivity is the same : 6 = 0.2. As expected, the phase diagram is no longer symmetric 
with dipolar hard-spheres : for a fixed total density p^, the coexisting phases are a phase 
rich in dipolar hard-spheres and a poor one, the later with a lower density, as shown by the 
positive slope of the tie lines. As an example, for p!j, = 0.51, 

one has p\ = 0.498 in the poor phase (x2 = 0.037) and pg = 0.523 in the rich one 
(xg = 0.96) . The second observation is that the lowest density of demixing (here nearly 
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the critical demixing point) decreases from 0.42 to ~ 0.395 with the dipoles. Both 
effects are rather moderate here since we took a low fi* for not departing too much from the 
conditions in which the population inversion occurs in the pure hard-spheres mixture 6|]. In 
this respect the dominant mechanism is here the non-additivity. In comparison, the demixing 
of a mixture of neutral and dipolar hard spheres with additive diameters should occur at 
much higher density [3], [l9|. For example the demixing of a mixture with X2 = 0.2 



occurs 



17l | at = 0.8 for p* = 1.8 when we have ~ 0.418 already for fi* = 1. This 



seems also consistent with the estimation made in ref. [16| from integral equations. Direct 
comparison with these studies is further complicated by the fact that they used different p*, 
the pressure instead of the density etc. Note also that this is not a pure demixing since the 
density is not constant {sS ]. 

With the same parameters of the model, we now consider the mixture confined in a slab. 



B. Fluid confined in a slab of width H = 3a 

1. Phase Diagram in confinement 

Fig. [2] shows the influence of the confinement on the demixing line. The noteworthy 
points are: 

- as in the bulk, the presence of dipoles favors the demixing (the lowest demixing density 
shifts from ^ 0.65 to p\ ^ 0.595). 

- as with non-additive hard spheres, the confinement stabilizes the mixture: the lowest 
demixing density shifts from p\ ^ 0.395 in the bulk to p^ ~ 0.595 for H = 3a. This effect 
depends of course on the pore width. Due to the non-additivity of the diameters between the 
particles and the pure hard core interaction between the latter and the walls, more particles 
can be packed in the pore than in a bulk slab of same width. This is consistent with the 
interpretation of the population inversion by an increased free volume in the pore for the 



minority species [10[. For a given width, the correspondence between a bulk state point and 
one in the pore is established by the equality of the chemical potentials in the bulk-open 
pore equilibrium, imposed through the GC simulation. This allows to discuss the position 
of a state point with respect to the corresponding demixing lines. For example, we found 
that a bulk state point (p!^ = 0.51,^2 = 0.02) corresponds to a single phase in equihbrium 



10 



with a single phase also in the pore with (p = 0.638, X2 = 0.06). Since the bulk mixtures 
with p\ = 0.51 demixes at X2 ~ 0.03 and the confined one with = 0.638 at X2 ~ 0.2, the 
distance of the state point (p = 0.638,^2 = 0.06) to the demixing boundary is larger than 
the corresponding one in the bulk. This shows that stabilization by confinement holds also 
for a mixture confined in an open pore. 



2. Connection with the density-induced population inversion 

Fig. [3] shows the population inversion in a pore in equilibrium with a bulk fiuid at fixed 
dipolar concentration = 0.02 and varying density. This figure corrects a slight inaccuracy 
in the similar plots in our previous papers after longer simulations the inversion 

is found to take place between = 0.53 and p\ = 0.54 and not right after p\ = 0.54. 
Comparison with the case of pure hard-spheres (figure 5 in ref. 6|) shows that the infiuence 
of the dipolar interaction is moderate. The starting point of the study of the field effect was 
indeed that the dipolar interaction should not disturb too much the density driven PINBI. In 



ref. jsl, we have shown that with pure non-additive hard-spheres , the population inversion 
occurs when the bulk fluid starts to become unstable (bulk state points slightly inside the 
demixed phase). We conflrm here this behavior when the non-additive mixture comprises 
dipolar species: as illustrated in Fig. HaJ where we show the path followed in the bulk fluid 
at the flxed composition = 0.02, the population inversion in the pore (points labeled 
2 and 3 in Fig. |4b] ) is indeed obtained for bulk state points (Fig. Ha|) slightly inside the 
bulk coexistence line. The state points in Fig. |4b] are however outside the pore demixing 
line. This conflrms that for this actual bulk composition, the population inversion is distinct 
from the phase separation in the pore. It corresponds to a sudden change in composition 
without coexisting phases: before and after the population inversion, there is only one phase 
in the pore but its composition jumps abruptly from nearly the bulk composition to the 
symmetric one 1 — This can be viewed as the equivalent in mixtures of the phenomenon 
of capillary condensation of a confined one-component fiuid. Note however that in order to 
observe the inversion as shown in Fig. [3l one needs to prepare the bulk fiuid in a metastable 
state, for these values of the physical parameters. 
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3. Effect of an external field 

Studying the effect of an external field should require mapping the {E* ,x^, p^) surface. 
Since this would require very lengthy simulations, we did that for selected values of the 
density. In Sec. IIII B 11 the value = 0.638 was considered to illustrate the effect of 
confinement on coexistence. It corresponds to the reference bulk state at which an important 
effect of the applied field on the composition of the confined fluid has been evidenced js, [gI. 
We thus show in Fig. [5] the effect of a uniform field (for two directions) on coexistence, the 
total density being fixed to p^ = 0.638. 

In a normal field {Ej_), the coexistence points (Fig. [Sa|) get closer at increasing field 
strength : starting from well separated zero-field values {p\ = 0.63, = 0.24) and (pg = 
0.645, ^2 = 0.75), the demixed domain shrinks as the field strength increases. This suggests 
a critical point beyond which the mixture becomes homogeneous. We also tentatively show 
f-PINBI points for the total density p^ = 0.51. They correspond to the bulk composition 
X2 = 0.02 (rightmost curve in figure 9 in 6|): {p^ = 0.51,^2 = 0.02) corresponds to (p = 
0.638, X2 = 0.06) in the bulk-pore equilibrium at zero-field. The coexistence curve E* — X2 
is for fixed at this value 0.638. It corresponds to a cut at the constant density p^ = 0.638 
of the {E* , X2, p!p) surface. Along the f-PINBI path now, the mean pore density p changes 
slightly from p ^ 0.62 to p ~ 0.638 for E* ^ — 3. It makes thus sense to show the f-PINBI 
data on the same plot as the E* — coexistence at fixed p^. We may then infer from their 
comparison that the field-induced population inversion occurs in a single-phase state in the 
pore for this specific bulk state. This scenario should still hold when the homogeneous bulk 
is at higher X2 since the tendency to demix is less pronounced in mixtures more symmetric 
in composition, a trend also favored by the mixing effect of the applied field. 

With the parallel field {E\\), the coexistence lines show the opposite trend as shown in 
Fig.l5b[ First, the demixed domain widens as E* increases. Second, the f-PINBI path crosses 
the demixing line. It is however difficult to draw a conclusion only from the position of the f- 
PINBI points since the variation of the mean pore density is now significant (p ~ 0.638 — 0.70 
for E* ^ — 1.5). From a study of the composition histograms in the GC simulation, we 
found that the two points closest to the coexistence line {E* = 1 and 1.25), the histograms 
are centered about a single value. At the intermediate value E* = 1.125, the histograms 
become bivariate. This shows that for = 0.02 there can exist a narrow range of field 
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strengths in which one has a demixed phase rather than a population inversion involving a 
single phase. 

In order to understand the physical mechanisms, it should be noted that the main effect 
of the applied field is to align the dipoles and decrease the potential energy in the pore by 
a contribution — 'Y^^i^ E. Its combination with the hard-core repulsions, the dipole-dipole 
interactions and the entropic effects minimizes the appropriate thermodynamic potential at 
equilibrium. The dipole-dipole interactions involve both attractive contributions between 
the dipoles that form a chain and repulsions between parallel chains. These positive con- 
tributions compete with the decrease in potential energy due to the field. The outcome of 
this combination is difficult to anticipate and depends on the kind of pore - closed or open - 
and on the field direction. Quite generally, the hard-core repulsion with the wall limits the 
length of the dipolar chains in the normal field and their number in the parallel one. 

In a closed pore (for which the E* — x\ coexistence curves at fixed are more directly 
relevant), the question is which situation - homogeneous or demixed mixture - leads to 
a minimum free energy. Starting from a demixed mixture (due to the non-additivity) at 
zero field, an increase of E±_ will favor mixing because only short dipolar chains can form. 
This also limits the decrease in entropy. The intra-chain attractive contributions that favor 
demixing are thus bounded and the interchain repulsion makes less efficient the effective 
repulsion between the neutral and the dipolar hard-spheres that drives the demixing. As a 
result a normal field favors the mixed state in the closed pore. In parallel field, formation of 
long chains parallel to the walls, is favored. An overall decrease of the potential energy can 
be anticipated due to these long chains. This can be directly checked for dipoles distributed 
on a regular lattice (see also fig. 10 in ref 30|). The fact that the field reinforces the 
trend towards demixing induced by the non-additivity shows that the dominant effect is the 
decrease in potential energy (this is consistent with the fact that the dipoles are naturally 
oriented parallel to the walls already at zero field). 

In the open pore, the grand potential is minimized also by a change in the number of 
particles (it varies during the GC simulation). Besides the mechanisms that act in the close 
pore, the filling of the pore by the dipoles is then energetically favorable. Accordingly, more 
and more dipoles enter the pore. With they form a homogeneous mixture with the 
hard-spheres, as for the closed pore. The population inversion is then progressive, with 
no interference with demixing. With £^11 the filling of the pore will interfere with field- 
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reinforced trend to demixing with a complex combination of both, including a change of the 
average density (no pure demixing). For example, in Fig. [5b] for the point corresponding to 
E* = 1.5, the mean number of dipoles N2 in the open pore is approximately 10% larger than 
the corresponding one at fixed density = 0.638. This larger number of dipoles that benefit 
from the decrease of potential energy due to the coupling with the field seems consistent 
with the preference for the population inversion rather than for the phase separation. For 
other parameters, the additional flexibility offered by the particle exchange with the bulk 



can lead to a great diversity of scenarios. 



as in the mesoscopic GC simulations [34 



'his would be especially true with a finite reservoir 



35|. 



C. Fluid in a slab of width H = 9a 

1. Phase Diagram in the pore 

Fig.[6]shows the phase diagram for a pore width H = 9a. The overall effect of confinement 
is the same as for H = 3a, but the situation is quantitatively intermediate between the latter 
case and the bulk. With the dipoles, the lowest density of demixing shifts from ^ 0.395 
in the bulk to ~ 0.435 in the slab. For further illustrating the lower confinement effect 
for H = 9a, we start with the same reference bulk state (p^ = 0.51, = 0.02) considered 
previously. Keeping the total total density in the pore fixed to p^ = 0.543 - that would be 
the average density (with X2 = 0.03) corresponding to the reference bulk state in the open 
pore-bulk equilibrium - and varying the dipole concentration Xg, the mixture demixes into 
a dipole poor phase with {p\ = 0.53, X2 = 0.07) and a rich one with (pg = 0.56, X2 = 0.95). 
At fixed total density p^ = 0.543, the mixture is homogeneous up to X2 ~ 0.05. Comparison 
with the discussion in Sec. IIIIB II shows that the stabilizing effect of confinement is less 
marked than for H = 3a. 



2. Effect of the external field 

Similarly to the case H = 3a, we start with the value of the total density p^ = 0.543 
considered above. The gross trends found for H = 3a are again observed (Fig. [7j) : a high 
normal field favors mixing while a parallel one favors demixing. The f-PINBI points are 
also shown, with the reservation about the representations discussed above. With E±, there 
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exists a range of field strengths {E* > 2) in which the population inversion goes clearly 
with a single phase (the PINBI point for E* = 2 which seems inside the demixed phase has 
an average density actually lower than 0.543). With E\\ the bulk-pore equilibrium leads to 
coexistence rather than a population inversion with a single phase (see also the discussion 
of the figure 14 in Gj). 

As the fluid with a fixed density = 0.543 is demixed at zero field almost over the 
entire composition range, the effect of the parallel field on the coexistence is difficult to see. 
We therefore considered a smaller density = 0.48 (Fig. |8]). At E* = 0, the corresponding 
dipolar concentration (in the poor phase), = 0.16, is comparable to that of the state 
point for which the field effect was studied for H = 3a. The shape of the coexistence lines 
are then similar to those for H = 3a (Fig. A qualitative difference is however that at 
low field strengths demixing is favored both for E± and E\\. But a sufficiently high normal 
field always favors mixing in sharp contrast to the effect of a parallel field. 



IV. CONCLUSION 



We may now summarize the main points of this study. Based on our previous work on a 
mixture of non-additive neutral and dipolar hard-spheres in the bulk and confined between 
hard walls, we have supplemented the results relative to the population inversion in an open 
pore with the coexistence curves for the bulk and the confined mixture. A combination of 
Monte Carlo simulations in the grand canonical ensemble for the pore-bulk equilibrium and 
in the Gibbs ensemble for the bulk-bulk and pore-pore equilibria made possible the use of the 
global density and composition of the bulk phase as the control variables. Accordingly, we 
have been able to discuss the interplay between population inversion and demixing with the 
help of these variables that are easily accessible experimentally. Our main motivation was to 
analyze further the connection between the phenomenon of population inversion in an open 
pore subject to an external field and the phase transitions in the bulk and in the confined 
fluid. Rather than an performing a generic study of the phase diagram, we considered a 
single value of the dipole moment, taken relatively low precisely for remaining close to the 
situation in which the phenomenon of density induced population inversion has been shown 
in previous studies. Other choices for the physical parameters or the study of transitions 
such as layering, segregation etc. might require different strategies than the GEMC/GC 
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combination considered here. 

Our main results concern the density driven population inversion and the field-induced 
one. At zero field, we observe that the presence of dipoles favors the demixing as in the 
bulk, but the confinement stabilizes the mixture. The density induced population inversion 
is then closely related to the bulk instability, as for the mixture of non-additive hard-spheres 
without dipoles. 

At non zero field, one may either consider its effect on phase separation at fixed total 
density, as in a closed pore, and both phase separation and population inversion, in an open 
one. The behaviors observed depend on the direction of the field. A normal field favors the 
mixing, and, for the open pore, a progressive population inversion without interference with 
demixing. In parallel field, demixing is favored. The population inversion interferes then 
with the phase separation. The result is a jump in the concentration of the dipolar species 
at the threshold field. Both effects are however difficult to separate, especially when the 
confinement is moderate. 

These observations can be rationalized largely in terms of the combination of the drop 
in electrostatic energy in the pore region for the dipolar species, the attractive intrachain 
dipolar interaction and the interchain repulsion, besides the entropy loss due to the alignment 
of the dipoles. Since these effects depend on the pore width and the direction of the field, 
their actual combination is difficult to anticipate. One should also consider the effect of 
the neutral hard-spheres whose overall effect is to decrease the electrostatic coupling. It is 
also difficult to draw general guidelines from the necessarily limited domain of the space 
of physical parameters than can be explored by simulation. Our previous studies of more 
complete versions of this model illustrate some aspects of such simulations {sgI . Is?! . Further 
studies of this model by semi-analytic methods should thus be useful at least for a gross 
estimation of the relevant range for these parameters (the same mixture without dipoles has 
recently been studied in the group by integral equations (38|) . 

From a practical point of view, we have shown that it is possible to select bulk states in 
which the bulk mixture is homogeneous and a field strength for which the confined mixture 
is also homogeneous but with an inverted population. This is precisely one of the impor- 
tant motivations of this study. In particular, the demixing can be avoided by a sufficient 
confinement as shown by the increase of the lowest demixing density when the pore width 
decreases. Finally, this analysis of the population inversion and the demixing phase tran- 
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sition has confirmed that a flexible effect of a uniform fleld can be observed by considering 
mixtures in the regime of high confinement. This can be paralleled with the case of bulk 
fluids or dielectric films in which a significant effect can be observed only when the field is 
inhomogeneous at a macroscopic scale [39| or the related studies of "electric bottles" [40 1. 
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Figures and caption 
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FIG. 1: Density-concentration demixing line of a bulk mixture of neutral and dipolar 
non-additive hard-spheres with 6 = 0.2. Full lines with crosses : /i* = 1, full lines with 
circles : fi* = (pure non-additive hard-spheres ). Tie lines are drawn between some pairs 

of coexistence points. 
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FIG. 2: Density-concentration demixing line of a mixture confined in a pore of width 
H = 3a. Meaning of symbols and lines as in Fig. [H 
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FIG. 3: Average densities in the confined mixture of neutral and dipolar non-additive 
liard-splieres as a function of the bulk density and fixed bulk concentration = 0.02. 
Other parameters as in Fig. [2] : = 3a, 6 = 0.2, /i* = 1; filled circles: adsorption of the 
dipolar hard-spheres, empty circles: desorption of the hard-spheres. The lines are guides to 

the eye. 




(a) Bulk path (demixing line shown in Fig. [1} (b) Corresponding path in the pore (demixing 

line shown in Fig. 

FIG. 4: Density-concentration trajectories in the phase diagram corresponding to the 
population inversion in Fig. |3l The bulk composition is = 0.02. The numbers label 
corresponding pore/bulk equilibrium state points. 
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(a) Normal field (b) Parallel field 

FIG. 5: Effect of tlie field strength on the demixing hnes at fixed total density = 0.638 
in the pore, corresponding for E* = to the bulk state (pt^ = 0.51, = 0.02). The 
associated f-PINBI path is shown by symbols. All the parameters are as in Fig. O 




FIG. 6: Phase diagram of a mixture of neutral and dipolar non-additive hard-spheres 
(/i* = 1) with 6 = 0.2 confined in a pore of width H = 9cr. 
Meaning of the lines and symbols as in Fig. [T] The point shown by a circle corresponds to 

the reference bulk state {p^ = 0.51, X2 = 0.02) 





(a) Normal field 



(b) Parallel field 



FIG. 7: Effect of the field strength on the coexistence line in a pore of width H = 9a 
fixed total density = 0.543. Meaning of the hnes and symbols as in Fig. |5l 



* 1 





0.2 0.4 0.6 0.8 1 



(a) Normal field 



(b) Parallel field 



FIG. 8: Same as Fig. [7] but at fixed total density = 0.48. 



